Measurement of radiative neutralino production 
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We perform the first experimental study with full detector simulation for the radiative 
production of neutralinos at the linear collider, at ^/s — 500 GeV and realistic beam 
polarizations. We consider all relevant backgrounds, like the Standard Model back- 
ground from radiative neutrino production. The longitudinal polarized beams enhance 
the signal and simultaneously reduce the background, such that statistical errors are 
significantly reduced. We find that the photon spectrum from the signal process can 
be well isolated. The neutralino mass and the cross section can be measured at a few 
per-cent level, with the largest systematic uncertainties from the measurement of the 
beam polarization and the beam energy spectrum. 

1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is a promising extension of the 
Standard Model of particle physics (SM) [T]. At a future Linear Collider (LC), the masses, 
decay widths, couplings, and spins of the new SUSY particles can be measured with high 
precision [2j. In particular, the lightest electroweak states like pairs of neutralinos, charginos, 
and slcptons, can be studied in the initial stage of the LC, with a center-of-mass energy 
y/s = 500 GeV, and a luminosity of order £ = 500 fb~^. The lightest visible SUSY state is 
a pair of radiatively produced neutralinos [3H6] 

e+ + e'^X? + X?+7- (1) 

The signal is a single high energetic photon, radiated off the incoming beams or off the ex- 
changed selectrons, and missing energy from the neutralinos. The main irreducible Standard 
Model background is photons from radiatively produced neutrinos e~^e~ — >■ vP"f^ see Fig.O 
Indeed, due to this large background, radiative neutralino production cannot be observed 
at LEP [5], not even for light or even masslcss neutralinos [7]. 

The discovery potential for the radiative production of light neutralinos is much better 
for the LC, since it provides high luminosity and the option of polarized beams [5]. In 
particular, it was shown that the electron and positron beam polarizations can be used to 
significantly enhance the signal and suppress the background from radiative neutrino pro- 
duction simultaneously I31[5J|B]. This applies in particular for scenarios where the lightest 
neutralino is mainly bino, i.e., has an enhanced coupling to the right sleptons. Right slepton 
exchange can be enhanced with positive electron and negative positron beam polarizations. 
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Figure 1: Contour lines of the tree- level cross section [SJIS] for radiative neutralino produc- 
tion, e'^e" XiXi7 (left) and radiative neutrino production e'^e^ — > vvj (right), at the 
LC for = 500 GeV and £ = 500 fb-\ for our benchmark scenario, Eq. ([2]). We applied 
regularization cuts on the photon angle | cos0^| < 0.99, and energy 5 < ^^-^[GeV] < 212. 



Since such an adjustment reduces the couplings to the W bosons, radiative neutrino produc- 
tion can be severely suppressed at the same time. In Fig. [1] we show the beam polarization 
dependence of signal and background, for our benchmark scenario with the relevant low 
energy SUSY parameters 

Ml = 103 GeV, Ah = 193 GeV, fi = 396 GeV, tan;3 = 10, mg^,^, = 125(190) GeV, (2) 

where the lightest neutralino m^o = 98 GeV has a bino component of Nu = 0.986. Already 
for a realistic beam polarization of {Pe- , Pe+) = (0.8, —0.3), we obtain large signal to back- 
ground ratios. The electroweak parameters of our scenario correspond to the widely studied 
SPSla' point [2]. Note that although the first generation squarks and gluinos below about 

1 TeV are excluded by recent LHC data [9], if their masses are roughly equal, the QCD 
sector of the SPSla' point is not relevant for our study. 

2 Experimental analysis 

Due to the large number of signal and irreducible background events in this channel, sys- 
tematic uncertainties of the detector measurements and of the beam parameters have to 
be considered. Therefore we perform an experimental study [TUHH], to demonstrate the 
potential of a LC for measuring the production cross section and the Xi mass. The analysis 
is made in the framework of the proposed International Linear Collider (ILC) [13] in full 
simulation of the International Large Detector concept (ILD) [14] , and is valid for the nom- 
inal ILC parameter set as given in the ILC Reference Design Report (RDR) [13]. The study 
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Figure 2: Photon energy distributions for e'^e'^ — > x5xi7 (black, solid), e'^e'^ — > vv^ (red 
dot-dashed), and e+e" ^ VD*-f (blue, dotted), at the LC for = 500 GeV, £ = 500 fb"\ 
and (Pg- , Pe+) = (0-8, —0.3), for scenario Eq. ([2|), with regularization cut | cos0^| < 0.99 [5]. 

which we shortly outline here [llj , is a model specific application of the model independent 
study to search for WIMP Dark Matter at the ILC P^[T5] . for further details see Ref. [TU]. 

2.1 Background and event generation 

Apart from the dominant irreducible SM background process of radiative neutrino pair 
production e+e" — )■ vv^{N)^ with up to three final state photons [N ~ 0, 1, 2), additional 
SM processes have been considered: Radiative Bhabha scattering e+e" — ?> e"'"e~7 with 
forward peaked final state electrons, leaving the detector through the beam pipe or hitting 
the forward calorimeters, as well as multi-photon final states e^e'' — >■ ^^(N)^, where only 
one of the emitted photons is properly reconstructed. 

For the event generation, Whizard [16l[l7] has been initialized with the ILC baseline 
parameter set for a -^/s = 500 GeV machine [T3|. The production of initial state radia- 
tion (ISR) in the leading logarithm approximation |16| is switched on, resulting in up to 
two predominantly soft and collinear additional photons. The beam energy spectrum used 
in the event generation has been calculated using GuineaPig [18] . 

2.2 Detector simulation 

For the detector simulation the Geant4 [19] based Mokka [20] simulation software has been 
used. The event reconstruction was performed with the PandoraPFA [21] reconstruction 
algorithm used within the MarlinReco [52] reconstruction framework for linear colliders. 
After simulation, two corrections are performed on the reconstructed event samples. In 
a first step, split uncharged electromagnetic clusters are iteratively merged with a cone 
based method to form higher level photon candidates. A second correction re-calibrates the 
detected photon energies with respect to energy lost in the cracks and inter-module gaps for 
cabling and readout of the ILD [TO] . 
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\Pe+\) = (0.8; 0.3) 



(|Pe-|; |i"e+l) = (0.8; 0.6) 




3.89 ± 0.8 (0.5) 

0.00 ± 1.1 (0.7) 

0.00 ± 1.2 (1.0) 

0.11 ± 1.3 (0.8) 



3.89 ± 0.3 (0.2) 

0.00 ± 0.8 (0.5) 

0.00 ± 0.8 (0.5) 

0.11 ± 1.0 (0.5) 



Table 1: Measured fully polarized cross sections a^R^Lj/ao ± (5^ [fb], see Eq. (O, for two 
different positron polarizations for scenario Eq. ([2|). The error is the squared sum of 
the statistical and the systematic errors. The values in brackets correspond to an improved 
polarization measurement of SP = 0.1%, instead of SP = 0.25%. 

2.3 Event sampling and selection cuts 

In order to enable an efficient study of many different models, the signal events were gen- 
erated by reweighting the fully simulated SM background samples. To avoid statistical 
correlation, the background samples have been divided into three statistically independent 
subsaniples. The first two subsamples are used for the background and signal contribution 
to the data. The third sample is parametrized to yield the predicted photon spectra. 

Since the lyvj SM background is indistinguishable on an event-by-event basis from the 
neutralino signal, the selection is tuned to isolate this SM background. An event is considered 
signal-like when it contains at least one photon with 



This signal definition ensures that the detected photons are within the tracking acceptance 
of the ILD detector to distinguish them from charged particles. The cut on the photon 
energy reduces the contributions from soft ISR, and excludes the massless neutrino final 
states from the radiative Z-return at photon energies of 241 GeV, see Fig. [5] 

In the further event selection, the following constraints are set to reject the dominant 
reducible SM backgrounds. To exclude hadronic and leptonic final states, a cut Evis — 

< 20 GeV is applied on the maximal exclusive energy, i.e., the total detected (visible) 
event energy excluding the selected photon. The maximal transverse track momentum is 
constrained to pT < 3 GeV. Low px tracks have to be allowed because of track overlays from 
e~^e~ pairs from the beamstrahlung background, and from multi-peripheral 77 — >■ hadrons 
events. Although the angular distribution of the final state fermions is strongly peaked in 
the forward direction, a strong activity in the sensitive volume is expected, due to the large 
cross section in the order of 5 x 10^ fb. The impact of both contributions has been studied for 
the ILD detector using Whizard for event generation [13]. The beamstrahlung spectrum 
has been simulated with Guinea Pig. On average 0.7 tracks from 77 processes and 1.5 
tracks from beamstrahlung background are expected per bunch crossing. Tight selection 
criteria on these tracks would therefore reduce the signal statistic strongly. Finally, the 
large Bhabha background is reduced by an identification of high energy electrons in the 
forward beam calorimctry. The selection efficiency of the viyj background is above 80% on 
average, while the multi-photon and Bhabha background is reduced to < 1%. The signal 
selection efficiency is about 92%. 



10 GeV <E^ < 220 GeV 



I cos 6*^,1 < 0.98. 



(3) 
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2.4 Cross section measurement and coupling structure 

The signal cross section is determined by a subtraction of the expected number of background 
events < Nb > from the number of observed data events No, 

(u u \ Nd-<Nb> 
aiP^-,P,+ ) = , (4) 

with the experimental luminosity C and the signal efficiency e. With four different polariza- 
tion configurations, the fully polarized cross sections cr{L.fl}, and hence the helicity structure 
of the coupling to the beam electrons, can be determined from [8] 
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(1 + P,-){1 - P,+ )aRL + (1 - Pe-)(1 + ^c+Ki?^ • (5) 



In Tab.[TJ we summarize the reconstructed cross sections and their corresponding errors. The 
luminosity of £ = 500 fb~^ is distributed to the odd (even) sign polarization configurations 
with C ~ 200 (50) fb~^ each. With an assumed systematic polarization measurement error 
of 5P = 0.25% [53], the cross section (Tjil can be determined to a level of about 20 (7)%, for 
an positron polarization of \Pe+ \ — 0.3 (0.6). Since the systematic uncertainty is dominated 
by the precision of the polarization measurement bP , the systematic error can be reduced 
by about a third if bP ~ 0.1% could be achieved, see Tab. [TJ At the ILC, the long term 
average of the luminosity weighted beam polarization can be determined from collision data 
to 0.17% for both electrons and positrons with |Pe+| — 0.6, dominated by the uncertainties 
from the polarimeters |24j . Under the assumption of uncorrelated polarization errors, a 
combination of the measured fully polarised cross sections cr{^ /^ j^i, yields an uncertainty of 
9.3 (6.0)% on the cross section for unpolarised beams for bP = 0.25% and \Pe+ \ = 0.3 (0.6), 
again systematically dominated by the polarimeters. 



2.5 Neutralino mass measurement 

The mass of the neutralino candidate is determined from a of the full data energy 

spectrum to template spectra with different masses. Depending on the degree of positron 
polarization, statistical precisions of 1.7 to 2.6 GeV can be obtained for an integrated lu- 
minosity of 500 fb~^. The total uncertainty on the candidate mass is dominated by the 
uncertainty in the beam energy spectrum, contributing with an error of 2.2 GcV to the 
mass determination, see Tab. [2l The influence of the beam energy spectrum, which dis- 
torts the energy spectrum of the signal photon, has been estimated by generation of signal 
spectra in a generic WIMP model [15^ for two different sets of beam parameters, RDR [13] 
and SB-2009 |25| . This is a conservative estimate, since the beam energy spectrum will be 
known to a higher degree than the difference between the sets. 
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±si&i. ±sys.{5E±5C) (total) [GcV] {P^-; P^+) 



97.7 
97.7 
97.7 



±2.65 ± 0.09 ± 2.20 (3.44) (0.8; 0.0) 

±2.07 ± 0.09 ± 2.20 (3.02) (0.8; -0.3) 
±1.70 ± 0.09 ± 2.20 (2.79) (0.8; -0.6) 



Table 2: Neutralino mass determined from a template comparison for an integrated lu- 
minosity of £ = 500 fb^^ and different beam polarizations. The systematic uncertainties 
comprise of the beam energy scale calibration {6E) and the beam energy spectrum {5C). 
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